Anharmonic IR spectra of solvated ammonium and aminium ions: resemblance between water and bisulfate solvations.

In this work, we analyze the vibrational spectra of ammonium, methylammonium, and dimethylammonium ions solvated by either water molecules or bisulfate anions using anharmonic vibrational algorithms. Rich and complicated spectral features in the 2700-3200 cm-1 region of the experimental spectra of these clusters are attributed to originate from strong Fermi resonance between hydrogen-bonded NH stretching fundamentals and NH bending overtones. Additional weaker bands around 2500-2600 cm-1 in solvated aminium ions are assigned to the combination tones involving the CH-NH (methyl-amino) rocking modes. Furthermore, the qualitative resemblance in band positions and spectral patterns between two-water-solvated and two-bisulfate-solvated cations suggest a common vibrational coupling scheme beneath the two seemingly different micro-solvation environments.

An ab initio anharmonic approach to IR, Raman and SFG spectra of the solvated methylammonium ion.

The methylammonium ion (CH3NH3+, or noted as MA-H+) is one of the smallest organic ammonium ions that play important roles in organic-inorganic halide perovskites. Despite the simple structure, the vibrational spectra of MA-H+ exhibit complicated features in the 3 µm region which are sensitive to the solvation environment. In the present work, we have applied the ab initio anharmonic algorithm at the CCSD/aug-cc-pVDZ level to simulate the IR and Raman spectra of the solvated methylammonium ion, MA-H+⋯X3, where X denotes the solvent molecules, to understand the Fermi resonance mechanism in which the overtones of NH bending modes are coupled with the fundamentals of NH stretching modes. The spectral features of the solvated clusters with proper solvent species resemble those observed in the perovskite crystal, indicating that they have similar solvation environments and hydrogen bond interactions. Therefore, a linkage between the gas-phase cluster models and the condensed-phase materials can be established, and our simulations and anharmonic analyses help in interpreting the spectral assignments of the observed IR and Raman spectra of perovskites reliably. Furthermore, we have extended this approach to the SFG spectra to demonstrate the selective appearance of bands depending on both the beam polarization configurations and the symmetry of vibrational modes.

Anharmonic Coupling Revealed by the Vibrational Spectra of Solvated Protonated Methanol: Fermi Resonance, Combination Bands, and Isotope Effect.

Intriguing vibrational features of solvated protonated methanol between 2400-3800 cm-1 are recorded by infrared predissociation spectroscopy. Positions of absorption bands corresponding to OH stretching modes are sensitive to changes in solvation environments, thus leading to changes in these vibrational features. Two anharmonic coupling mechanisms, Fermi resonance (FR) contributed by bending overtones and combination band (CB) associated with intermolecular stretching modes, are known to lead to band splitting of OH stretching fundamentals in solvated hydronium and ammonium. Theoretical analyses based on the ab initio anharmonic algorithm not only well reproduce the experimentally observed features but also elucidate the magnitudes of such couplings and the resulting interplay between these two mechanisms, which provide convincing assignments of the spectral patterns. Moreover, while the hydroxyl group plays the leading role in all the above-mentioned features, the role of the methyl group is also analyzed. Through the H/D isotope substitution, we identify overtones of the methyl-hydroxyl rocking modes and their participation in FR.

Strong Fermi Resonance Associated with Proton Motions Revealed by Vibrational Spectra of Asymmetric Proton-Bound Dimers.

Infrared spectra for a series of asymmetric proton-bound dimers with protonated trimethylamine (TMA-H+ ) as the proton donor were recorded and analyzed. The frequency of the N-H+ stretching mode is expected to red shift as the proton affinity of proton acceptors increases. The observed band, however, shows a peculiar splitting of approximately 300 cm-1 with the intensity shifting pattern resembling a two-level system. Theoretical investigation reveals that the observed band splitting and its extraordinarily large gap of around 300 cm-1 is a result of strong coupling between the fundamental of the proton stretching mode and overtone states of the two proton bending modes, that is commonly known as Fermi resonance (FR). We also provide a general theoretical model to link the strong FR coupling to the quasi-two-level system. Since the model does not depend on the molecular specification of TMA-H+ , the strong coupling we observed is an intrinsic property associated with proton motions.

Vibrational spectroscopic signatures of hydrogen bond induced NH stretch-bend Fermi-resonance in amines: The methylamine clusters and other N-H⋯N hydrogen-bonded complexes.

The appearance of multiple bands in the N-H stretching region of the infrared spectra of the neutral methylamine dimer and trimer is a sign of NH bend-stretch anharmonic coupling. Ab initio anharmonic calculations were carried out in a step-wise manner to reveal the origin of various bands observed in the spectrum of the methylamine dimer. A seven-dimensional potential energy surface involving symmetric and asymmetric stretching and bending vibrations of both the hydrogen bond donor and the acceptor along intermolecular-translational modes was constructed using the discrete variable representation approach. The resulting spectrum of the dimer shows five bands that can be attributed to the symmetric stretching (νsym D), asymmetric stretchin (νasym D), and bending overtone (2νbend D) of the donor moiety. These appear along with the combination band arising out of bending vibrations of the donor and acceptor (νbend D + νbend A) and with the combination of the intermolecular translational mode over the donor bending overtone (νtrans + 2νbend D). The spectrum of the trimer essentially consists of all the features seen in the dimer with marginal changes in band positions. The analysis of the experimental spectra based on the two-state deperturbation model and ab initio anharmonic calculations yield a matrix element of about 40 cm-1 for the N-H bend-stretch Fermi resonance coupling. In general, the IR spectra of the hydrogen-bonded amino group depict three sets of bands that arise due to bend-stretch Fermi resonance coupling.

Anharmonic coupling behind vibrational spectra of solvated ammonium: lighting up overtone states by Fermi resonance through tuning solvation environments.

Studies on the vibrational spectra of various ammonium-centered clusters under different solvation environments have raised interest over the last thirty years. The gas-phase infrared photodissociation spectroscopy (IRPD) experiments showed that these NH4+Xn clusters exhibit rich spectral features from 2600 to 3400 cm-1. In this work, we have simulated the vibrational spectra and analyzed couplings among vibrational quantum states in the aforementioned frequency range using ab initio anharmonic algorithms. Originating from the anharmonic couplings between NH stretching fundamentals and bending overtones, Fermi resonance (FR) is a common feature in these spectra, and its extent is determined by the magnitude of couplings and the energy matching conditions between relevant states, which are governed by the proton affinity, number, and bonding configuration of the solvation species. For weakly bound clusters consisting of rare gas atoms, FR is insignificant but not negligible; for strongly bound clusters, such as ammonium-water clusters, the hydrogen-bonded NH stretching fundamentals redshift and reach a better resonance condition, and thus light up the bending overtones as prominent FR bands. Our simulated spectra are in good agreement with previous experimental reports of these ammonium-centered clusters and provide a better understanding of the vibrational coupling behind the spectra of the NH stretching region.

Vibrational spectroscopy of protonated amine-water clusters: tuning Fermi resonance and lighting up dark states.

Strong coupling between stretching fundamentals and bending overtones of vibrational modes, known as Fermi resonance (FR), has been observed for proton motions in the protonated trimethylamine-water cluster. To investigate the role of FR, we examined the vibrational spectra of other three protonated ammonia/amine-water clusters, including the NH4+ ion and its mono- and di-methylated analogues, respectively, with and without argon tagging. In these systems, a simple frequency-scaled harmonic oscillator model will predict only one strong band between 2600 and 3200 cm-1 uniquely due to the hydrogen-bonded NH stretching fundamental for a given conformer. In the experimental vibrational spectra, however, multiple sharp bands were observed. Such a discrepancy often leads to the notions of the coexistence of multiple conformers and/or the appearance of an overtone state as a result of FR. In this work, we applied a discrete variable representation (DVR) implementation of ab initio anharmonic algorithms and demonstrated how one N-H+ stretching fundamental can lead to multiple bands as a result of intrinsic anharmonic couplings. A prominent effect of tuning these FR bands and lighting up dark overtone states in this wide frequency range was investigated by changing the number of methyl groups in the protonated amine moiety. The effect of Ar-tagging was also analyzed and decent agreement between the experimental and simulated spectra certified the above-mentioned simple pictures. We also found that the coupling constant for trimethylamine is the largest among these protonated amine-water clusters, and the overall coupling strength decreases as the hydrogen-bonded NH stretching frequency redshifts in the order of dimethylamine, methylamine, and ammonia.

Improved agreement between experimental and computational results for collision-induced dissociation mass spectrometry of cation-tagged hexoses.

To model the collision-induced dissociation mass spectrometry (CID-MS) of Na+-tagged hexoses, it is not only required to perform an extensive sampling of the conformational space as addressed in our previous work [Huynh et al., Phys. Chem. Chem. Phys., 2018, 20(29), 19614-19624], it is also necessary to apply a sufficiently reliable quantum chemical method to describe the relative energetics of the reactions. In this work, the ring-opening via hemiacetal scission and the competing dehydration pathways have been re-evaluated at the MP2 level. The results show that previous studies at the B3LYP level display a systematic underestimation of the dehydration barriers by about 40 kJ mol-1 on average while the ring-opening barriers are reasonably described. We further illustrate that, in the present case, it is not enough to only look at the energetics: although MP2 results indicate that the ring-opening pathways of the considered hexoses have lower barriers than the dehydration pathways, the contributions of the partition functions to the rate constants render the dehydration for some structures to be more favorable at elevated temperatures. Via benchmark calculations against single-point CCSD(T) results at the example of small organic model molecules, we demonstrate that the MP2 data for the studied reactions are, alongside calculations at the M06-2X level, trustworthy. In addition, we have analyzed for the same set of model molecules how the fraction of exact exchange in the applied exchange-correlation functional affects the reaction barriers and the charge distributions. While MP2 (and M06-2X) calculations indicate an increased charge localization when going from the initial state to the dehydration transition state, this phenomenon is not seen in the B3LYP calculations, which likely is the origin of the discrepancy between the B3LYP and MP2 dehydration barriers. As the variation of the charge localization during a hemiacetal scission reaction is comparably small, B3LYP performs sufficiently well for such reactions.

First-principles study on sum-frequency generation spectroscopy of methanol adsorbed on TiO2(110) surface: Effects of substrate and molecular coverages.

In this work, starting from the general theory of sum-frequency generation (SFG), we proposed a computational strategy utilizing density functional theory with periodic boundary conditions to simulate the vibrational SFG of molecules/solid surface adsorption system. The method has been applied to the CH3OH/TiO2(110) system successfully. Compared with the isolated molecule model, our theoretical calculations showed that the TiO2 substrate can significantly alter the second-order susceptibilities of a methanol molecule which is directly related to the SFG intensity. In addition, the SFG spectra have obvious changes while the methanol coverage increases, especially for the OH vibration peaks. Our theoretical spectra agree reasonably well with experimental measurements at 1 ML coverage, and an interesting peak which is absent in the theoretical spectra is tentatively assigned to some CH3 stretch vibration of methanol adsorbed on the oxygen vacancy of TiO2.

The dp type π-bond and chiral charge density waves in 1T-TiSe2.

Based on the atomic electronic configuration and Ti-Se coordination, a valence bond model for the layered transition metal dichalcogenide (TMDC) 1T-TiSe2 is proposed. 1T-TiSe2 is viewed as being composed of edge-sharing TiSe4-plaquettes as TiSe2-ribbon chains in each layer via a directional valence shell electron distribution as chemical bonds, in contrast to the conventional layer view of face-sharing TiSe6-octahedra. The four valence electrons per Ti in the hybridized dsp2-orbitals of square coordination form σ-bonds with the four nearest neighbor Se atoms in the chain. The electrons in the lone pair of the Se-4pz orbital are proposed to form a dp type π-bond via side-to-side orbital overlap with the empty Ti-3dxz/3dyz orbitals within each chain, which is positively supported by quantum chemistry calculations. A study of electron energy loss spectroscopy (EELS) with transmission electron microscopy (TEM) for 1T-TiSe2 is presented to show an energy loss near ∼7 and ∼20 eV, which confirms the existence of collective plasmon oscillations with the predicted effective electron numbers for the π- and (π + σ)-bond electrons, respectively.

Theoretical study of nitrogen-doped graphene nanoflakes: Stability and spectroscopy depending on dopant types and flake sizes.

As nitrogen-doped graphene has been widely applied in optoelectronic devices and catalytic reactions, in this work we have investigated where the nitrogen atoms tend to reside in the material and how they affect the electron density and spectroscopic properties from a theoretical point of view. DFT calculations on N-doped hexagonal and rectangular graphene nanoflakes (GNFs) showed that nitrogen atoms locating on zigzag edges are obviously more stable than those on armchair edges or inside flakes, and interestingly, the N-hydrogenated pyridine moiety could be preferable to pure pyridine moiety in large models. The UV-vis absorption spectra of these nitrogen-doped GNFs display strong dependence on flake sizes, where the larger flakes have their major peaks in lower energy ranges. Moreover, the spectra exhibit different connections to various dopant types and positions: the graphitic-type dopant species present large variety in absorption profiles, while the pyridinic-type ones show extraordinary uniform stability and spectra independent of dopant positions/numbers and hence are hardly distinguishable from each other. © 2018 Wiley Periodicals, Inc.

Pyridine vs N-Hydrogenated Pyridine Moieties: Theoretical Study of Stability and Spectroscopy of Nitrogen-Contained Heterocyclic Aromatic Compounds and Graphene Nanoflakes.

Nitrogen is one of the most common heteroatom appearing in heterocyclic aromatic compounds (HACs) as well as the frequently applied dopant in graphene nanoflakes/nanoribbons. The pyridine moiety is an intuitive and stable common feature of these compounds; but interestingly, using density functional theory calculations, we found that the N-hydrogenated pyridine moiety could be even more stable in large HACs and in N-doped graphene nanoflakes considering their formation reaction energies. The hydrogenation reaction of the pyridine moiety was calculated to be exothermic for models of four and more fused aromatic rings with specific substitutional positions of nitrogen. This theoretical investigation provides energetic and spectroscopic hints to the existence of the N-hydrogenated pyridine moiety under proper conditions.

Generalizability Theory With One-Facet Nonadditive Models.

In generalizability theory (G theory), one-facet models are specified to be additive, which is equivalent to the assumption that subject-by-facet interaction effects are absent. In this article, the authors first derive estimators of variance components (VCs) for nonadditive models and show that, in some cases, they are different from their counterparts in additive models. The authors then demonstrate and later confirm with a simulation study that when the subject-by-facet interaction exists, but the additive-model formulas are used, the VC of subjects is underestimated. Consequently, generalizability coefficients are also underestimated. Thus, depending on the nature of interaction effects, an appropriate model, either additive or nonadditive, should be used in applications of G theory. The nonadditive G theory developed in this article generalizes current G theory and uses data at hand to determine when additive or nonadditive models should be used to estimate VCs. Finally, the implications of the findings are discussed in light of an analysis of real data.

The role of the nπ* 1A(u) state in the photoabsorption and relaxation of pyrazine.

The geometric, energetic, and spectroscopic properties of the ground state and the lowest four singlet excited states of pyrazine have been studied by using DFT/TD-DFT, CASSCF, CASPT2, and related quantum chemical calculations. The second singlet nπ* state, (1)A(u), which is conventionally regarded dark due to the dipole-forbidden (1)A(u)←(1)A(g) transition, has been investigated in detail. Our new simulation has shown that the state could be visible in the absorption spectrum by intensity borrowing from neighboring nπ* (1)B(3u) and ππ* (1)B(2u) states through vibronic coupling. The scans on potential-energy surfaces further indicated that the (1)A(u) state intersects with the (1)B(2u) states near the equilibrium of the latter, thus implying its participation in the ultrafast relaxation process.

The visible spectrum of zirconium dioxide, ZrO2.

The electronic spectrum of a cold molecular beam of zirconium dioxide, ZrO(2), has been investigated using laser induced fluorescence (LIF) in the region from 17,000 cm(-1) to 18,800 cm(-1) and by mass-resolved resonance enhanced multi-photon ionization (REMPI) spectroscopy from 17,000 cm(-1)-21,000 cm(-1). The LIF and REMPI spectra are assigned to progressions in the Ã(1)B(2)(ν(1), ν(2), ν(3)) ← X̃(1)A(1)(0, 0, 0) transitions. Dispersed fluorescence from 13 bands was recorded and analyzed to produce harmonic vibrational parameters for the X̃(1)A(1) state of ω(1) = 898(1) cm(-1), ω(2) = 287(2) cm(-1), and ω(3) = 808(3) cm(-1). The observed transition frequencies of 45 bands in the LIF and REMPI spectra produce origin and harmonic vibrational parameters for the Ã(1)B(2) state of T(e) = 16,307(8) cm(-1), ω(1) = 819(3) cm(-1), ω(2) = 149(3) cm(-1), and ω(3) = 518(4) cm(-1). The spectra were modeled using a normal coordinate analysis and Franck-Condon factor predictions. The structures, harmonic vibrational frequencies, and the potential energies as a function of bending angle for the Ã(1)B(2) and X̃(1)A(1) states are predicted using time-dependent density functional theory, complete active space self-consistent field, and related first-principle calculations. A comparison with isovalent TiO(2) is made.

Symmetry forbidden vibronic spectra and internal conversion in benzene.

The spectra of symmetry-forbidden transitions and internal conversion were investigated in the present work. Temperature dependence was taken into account for the spectra simulation. The vibronic coupling, essential in the two processes, was calculated based on the Herzberg-Teller theory within the Born-Oppenheimer approximation. The approach was employed for the symmetry-forbidden absorption/fluorescence, and internal conversion between 1(1)A(1g) and 1(1)B(2u) states in benzene. Vibrational frequencies, normal coordinates, electronic transition dipole moments, and non-adiabatic coupling matrix elements were obtained by ab initio quantum chemical methods. The main peaks, along with the weak peaks, were in good agreement with the observed ones. The rate constant of the 1(1)A(1g)← 1(1)B(2u) internal conversion was estimated within the order of 10(3) s(-1). This could be regarded as the lower limit (about 4.8 × 10(3) s(-1)) of the internal conversion. It is stressed that the distortion effect was taken into account both in the symmetry-forbidden absorption/fluorescence, and the rate constants of internal conversion in the present work. The distortion effects complicate the spectra and increase the rate constants of internal conversion.

A theoretical study on the spectroscopy and the radiative and non-radiative relaxation rate constants of the S(0)(1)A(1)-S(1)(1)A(2) vibronic transitions of formaldehyde.

We have carried out a close examination on the mathematical treatments and the first-principle computations concerning the vibronic transitions between the S(0)(1)A(1) and the S(1)(1)A(2) states of formaldehyde. The simulation of absorption spectrum was presented with peak intensities calculated according to vibronic-coupled transition dipole moments and Franck-Condon factors. The radiative and non-radiative transition rate constants from the excited to the ground states were calculated with formulas based on Fermi's golden rule. It is concluded that our simulated absorption spectrum between 300 and 360 nm, as well as the estimated relaxation rate constants, showed good agreements with experimental reports.

One- and two-photon absorption properties of diamond nitrogen-vacancy defect centers: A theoretical study.

The negatively charged nitrogen-vacancy defect center, (NV)(-), in diamond has been investigated theoretically for its one- and two-photon absorption properties involving the first excited state with the (3)A(2)-->(3)E transition. Time-dependent density functional theory (TD-DFT), configuration interaction with single excitation (CIS), and complete active space self-consistent field (CASSCF) were employed in this investigation along with the 6-31G(d) basis set. Diamond lattice models containing 24-104 carbon atoms were constructed to imitate the local environment of the defect center. TD-DFT calculations in large molecular cluster models (with 85 or more carbon atoms) predicted the vertical excitation energy quite consistent with the experimental absorption maximum. CASSCF calculations were feasible only for small cluster models (less than 50 carbon atoms) but yielded one-photon absorption (OPA) and two-photon absorption (TPA) cross sections somewhat larger than the experimental values obtained with linearly polarized incident light [T.-L. Wee et al., J. Phys. Chem. A 111, 9379 (2007)]. CIS calculations in larger cluster models showed a systematic overestimation of the excitation energy while just slightly underestimated the OPA cross section and overestimated the TPA cross section. The agreements between calculations and measurements suggest that the computational approaches established in this work are applicable to explore the optical properties of related defect centers in diamond as well.

Symmetric double-well potential model and its application to vibronic spectra: studies of inversion modes of ammonia and nitrogen-vacancy defect centers in diamond.

In this paper, we have studied the vibronic transitions between two symmetric double-well potentials by proposing a model Hamiltonian consisting of a harmonic oscillator and a parturition described by a Gaussian function that leads to a double minima potential with a barrier between the two energy minima. Making use of the contour integral form of Hermite polynomials, we present a new formula that can calculate Franck-Condon factors of the system rigorously. The simulated vibronic spectra of ammonia and the negatively charged nitrogen-vacancy center in diamond are presented as an application of the formula.

Calculation of the vibrationally non-relaxed photo-induced electron transfer rate constant in dye-sensitized solar cells.

In this paper we shall show how to calculate the single vibronic-level electron-transfer rate constant, which will be compared with the thermal averaged one. To apply the theoretical results to the dye-sensitized solar cells, we use a simple model to describe how we model the final state of the electron-transfer process. Numerical calculations will be performed to demonstrate the theoretical results.